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Abstract 

In this paper we show that a recent hydrodynamical model of superfluid turbulence 
describes vortex density waves and their effects on the speed of high-frequency second 
sound. In this frequency regime, the vortex dynamics is not purely diffusive, as for low 
frequencies, but exhibits ondulatory features, whose influence on the second sound is here 
explored. 

Hydrodynamics of superfluid turbulence, characterized by a tangle of quantized vortex 
lines [HI2], aims to describe the couplings between pressure and heat perturbations, and the 
vortex density dynamics [3]-[7j. Such hydrodynamics is a lively topic, with recent emphasis, 
for instance, on nonlinear features such as the influence of intense heat pulses on the vortex 
tangle [8] -[10], or in multi-scale formulations allowing to eliminate the fast processes to derive 
evolution equations for the slow processes In the present paper we will focus our attention 
on a different topic, namely, the behaviour of linear heat waves and vortex density waves in 
the high-frequency domain. This has not received much previous attention, but we think it 
is worthwhile because, whereas at low frequencies the behaviour of the vortex distortions is 
mainly diffusive, at high enough frequencies, the vortex lines behave as an elastic medium, 
and are able to propagate waves by themselves, i.e. they behave like a viscoelastic medium: 
diffusive at low frequencies, elastic at high frequencies. Their effect on second sound propaga- 
tion is of much interest to provide a suitable physical interpretation of the experimental data 
based on second sound. 

Collective vortex waves in rotating superfluids have been studied in depth for a long time 
|12j , and the so-called Tkachenko transverse elastic waves in the vortex arrays with crystalline 
order arising in rotating cylinders have been theoretically discussed since 1966. However, 
vortex density waves in counterflow situations have not been studied, up to our knowledge. 
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Usually, one considers homogeneous vortex tangles in counterflow situation (i.e. under 
vanishing barycentric speed) and the evolution of L, the vortex-line density, is assumed to be 
the well-known Vinen's equation [TJ EJ EJ [13] 

^ = a L = AqLW - Bh\ (1) 

where q is the absolute value of the heat flux q, A and B parameters linked to the dimensionless 
coefficients appearing in Vinen's equation by the relations A = a v /p s Ts and B = Kj3 v [7], 
p s being the density mass of the superfiuid component, T the absolute temperature, s the 
entropy, k = h/m the quantum of vorticity (m the mass of the 4 He atom and h Planck's 
constant, k ~ 9.97 10 -4 cm 2 /s). 

This equation is used to describe homogeneous turbulence. When the vortex line density 
is inhomogeneous, such an equation must be generalized incorporating nonlocal terms. In [TJ, 
a set of evolution equation for e, q and L, with e being the energy density, was proposed with 
special attention to their consistency with the second law of thermodynamics through the Liu's 
procedure. Such equations were [7] 

pe + V • q = 0, 

q + CoVr + xoVL = -KLq, (2) 
L + V • (u q) = o L . 

The coefficient K is related to the Hall-Vinen friction coefficient Buy as K = t^kBhv, u o 
and xo are suitable coefficients whose physical meaning we want to explore and Co is the ratio 
of the heat conductivity of the fluid over the relaxation time of the heat flux, which is given 
by (KL)^ 1 . A further term proportional to L 3 / 2 q may be added to the right-hand side of 
the equation (2b) for the heat flux, describing a dry friction term. Here, we will omit it for 
the sake of simplicity [6j [7J . For vanishing xo and v o , the two first equations describe second 
sound (temperature waves) with speed V" 2 2 = -^r, with cy the specific heat per unit mass at 
constant volume, and the third equation reduces to Vinen's equation (1). In the steady state, 
the second equation reduces to 

The coefficient of the first term is just the heat conductivity already mentioned, and the 
second term describes a coupling between the vortex line density and the heat flux, analogous 
in some way to the coupling between concentration gradient and heat flux in usual fluids (the 
so-called Soret effect). The term in VL could be related to the internal energy density of the 
tangle, which is given by eyL, with ey = (p s K 2 /4ir) \n.(c/ a^L 1 / 2 ) the energy per unit length 
of the vortices. In more specific terms, it may be related to the tangle contribution to the 
total pressure, as given in (10). Since in this paper we are concerned with high frequency and 
short wavelength perturbations, the contribution of VL will be relevant, in contrast with usual 
situations at long wavelengths. 

From (2) it is seen that when q is neglected and the system is isothermal, and in the linear 
approximation, one obtains from (2b) and (2c) a reaction-diffusion equation for L of the form 

L = ——VL + a L , (4) 



whose diffusion coefficient D is given by D = -j^-- 
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If, instead, q is not neglected in (2b) and still in the isothermal assumption - for the sake 
of simplicity - we obtain 



L + KLqL = v oX oV L + KL Q a L + & L . (5) 

At high frequencies, i.e. for uj S> KLq and much higher than the inverse of the characteristic 
time of the vortex destruction and formation as described by ctl, namely uj 3> 2BLq — ^AqL^ 2 , 
(5) becomes 

L « z^oXoV 2 L, (6) 

that is, we get vortex waves with speed 

vlo = voXo- (7) 

This outlines the physical features of the coefficients vqXo as their product is seen to be related 
both to the diffusion coefficient and to the speed of vortex waves. 

Elastic compressional waves in the vortex density obtained in (6) may also be interpreted 
in a more mechanistic way. First of all, we may interpret the term vqc{ in equation (2c) 
as z/oq = Lvl, with vl an effective velocity of the tangle because the divergence term in the 
balance equation (2c) may be interpreted as a flux of vortex lines moving with peculiar velocity 
vl- This effective velocity will have an evolution equation of the type 

Pzff^jf = ~VPv, (8) 

where p e ff is an effective density of the vortices and py is the vortex contribution to the 
pressure. In fact, vortices themselves do not have mass, but they have indirectly associated 
an inertia due to the mass of the surrounding rotating fluid. By following Sonin (page 94 of 
Ref. 12) we estimate this mass per unit volume as p s Lr 2 , with p s the density of the superfluid 
component, which is the one which participates in the rotation around the vortices of the fluid 
and r a characteristic radius of the zone in which the superfluid is affected by the motion of the 
vortex; this will be of the order of the average vortex separation, namely L -1 / 2 ; in this way, 
we have p e ff oc p s . This is the reason that p s is the density which appears in the coefficient 
ey appearing in the expression for the vortex density of the tangle, as given in the paragraph 
below equation (3). 

Then, we may combine equation (8) with equation (2c), which at high enough frequency, 
when the influence of the production-destruction term is negligible, yields 

— — = -L V • — = V p v . (9) 

dt 2 dt p eff 

The total pressure of the turbulent superfluid has the form (equation (4.16) of Ref. 7) 

P = p* + e v L, (10) 

p* being the pressure of the bulk superfluid and eyL the contribution of the tangle, with ey the 
energy per unit length of the vortices. Combination of (9) and (10) and taking into account 
that p e ff oc p s , yields 



d 2 L 
dt 2 



L ^)V 2 L. (11) 

Peff J 



3 



Then the velocity of the vortex waves will be 



-i=^o^oc^lnf^-V (12) 

The combination Lqk? is similar to the combination Ok appearing in the velocity of inertial 
waves in the vortex array in rotating superfluids [12], if one replaces Lq = 20,/ k, O being the 
angular speed of the rotating cylinder. 

We now go to the diffusion coefficient D corresponding to equation (4). Since we have 
obtained an estimation of foXo appearing in (6), we will take advantage of it to obtain an 
expression for D. We get 

D = ^oc^^ln(') (13) 
KL B hv 4tt V«o£ 1/2 y 

Indeed, it was known that D oc k on dimensional grounds and on some numerical sim- 
ulation [14] . Here, we have a more explicit expression, based on a more microscopic model. 
Unfortunately, our derivation cannot set precisely the proportionality constant in (13) due to 
the ambiguity in the mass density p e fj associated to the vortices (more precisely, to the su- 
perfluid associated to the vortices). The analysis of the high-frequency regime may be indeed 
rewarding in new results which are not available if the study is limited to the low-frequency 
regime. 

Up to here we have considered an isothermal situation. According to (2), the behavior of 
q and L under non-isothermal conditions is connected to the behavior of the field T in such 
a way that the complete study of the system (2) is required. In [7J two of us have already 
tempted to solve this more complicated situation, but now to the light of which we have said 
before, some more explicit conclusions may be done with a clearer physical meaning in the 
high frequency regime, which complements the information obtained at low frequencies. 

As in f?]. expressing the energy in terms of the vortex line L and the temperature T and 
linearizing the two contributions ctl and Lq around the stationary solutions, already found in 

0, 

q = qD = (gio,0,0), L = L = ^[qio} 2 , T = T (x) = T* - (14) 
the system (2) assumes the following form 
pcyT + peLL + V • q = 0, 

q + CoVT + xoVL = -A'[L q+q (L-L )], (15) 
L + vqV ■ q = 



2BL - UqioL'J 2 



Q /O 

(L - L ) + Aq w L ' (qi - q 10 ) 



Now, if we suppose the propagation of harmonic plane waves of the form 

T = T (x) + fe i ( Ra - x -^ 

q = q + faHKn-x-ut) ( 16 ) 
L = Lq + £e i ( j?n - x - w *), 

where K = k r + ik s is the complex wave number, oj is the frequency and n the unit vector in 
the direction of the wave propagation, then we have the propagation of waves (vortex waves 
and heat waves) which move at the same speed W2 = jr- given by the combination of the second 
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sound and vortex waves velocities. As in [7j, we are assuming that the quantities oversigned 
by a tilde are small ones, and whose product can be neglected. 

Inserting (16) in the linearized system (15) we obtain the following algebraic system for 
the quantities T, L and q 

-[pc v ] uf - [peL}ouL + • n = 0, 

(-w - iN x )q + K[( ] fn + (K[ Xo } n - iN 3 q ) L = 0, (17) 
(-W - iN 2 )L + K[u ] q_ • n + iN± qi = 0, 

where 

N 1 = KL Q , N 2 = 2BL -^ALl /2 qw , N 3 = Kq 10 , iV 4 = Aq 10 L 3 Q /2 , 

and the subscript 0, which will be deleted from now on, denotes quantities referring to the 
unperturbed states. 

Now, let consider that the direction of the wave propagation is collinear to the initial heat 
flux, i.e. n = (1,0,0), then the system (17) becomes 

—pcyuT — peiojL + Kq\ = 0, 
(_ w _ iJVi)& + K Cot + (K Xo - iN 3 ) L = 0, 

(-w - iNt)q 2 = 0, (18) 
l-u - iNi)q 3 = 0, 
. (-u - iN 2 )L + {Kvq + iN A ) qi = 0. 

In the hypothesis of high-frequency waves, which means u S> N%, u ^> N 2 and |ET[ S> 

(3/2 \ 
^x^~i A9l °ut° — ) ' P rev i° us algebraic system (18) becomes 

—pcyufF — pej^uL + k r qi = 0, 
-uqi + k r C f + KxoL = 0, 

— ujL + k r i>oqi = 0, (19) 
-ujq 2 = 0, 
-ujq 3 = 0. 

The above system has nontrivial solutions if and only its determinant is zero, which corre- 
sponds to the following dispersion relation 

w 2 2 = V 2 2 (l-v Q pe L ) + vl , (20) 

where w 2 = u/k r is the speed of the wave, V 2 = Co/pcv is the second sound speed in the 
absence of vortices and = Xo v o 1S t ne speed of the vortex wave, which we have found in 
(7). We recall that all three fields T, L and q vibrate with the same speed w 2 given by (20) 
in such a way that each field contributes to the vibrations of the other two. If we try to read 
the relation (20) in terms of the second sound, we note that the vortex vibrations modify this 
second sound speed through the two contributions —V 2 vqP€l and v^, the latter due to the 
presence of the vortex waves and the former due to the reciprocal existence of two waves. The 
same conclusion may be achieved reading the relation (20) in terms of the vortex waves. The 
correction for the speed of the second sound is not important, because V 2 is of the order of 20 
m/s near 1,7 K [HUH], whereas, for Lq = 10 6 cm -2 , and according to the estimation (12), the 
speed of vortex density waves would be of the order 0, 25 cm/s, much lower than V 2 . 

In the earlier analysis of the system (15), we have only considered the terms of the equations 
in which to and k r appear such as we have also neglected the term k s relative to the dissipation 
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of the wave. Now, we assume that the quantities N±, N2, N3 and N4 are coefficients small 
enough to assume them as perturbations of the physical system. This is reasonable at high- 
frequencies, since we have assumed that u S> Nx, u> N2 and \K\ 3> max f , Aqw J^° — j ■ 

Thus, the contributions of these coefficients can modify the speed of the wave in a small 
quantity 5 and the imaginary part of the wave number, k s . 

Therefore, let assume that the speed of the wave has the following expression 

w = = w 2 + 5, (21) 

for which substituting it in the dispersion relation, i.e. the equation obtained imposing that 
the determinant relative to the system (15) vanishes, we obtain at the lower order the relation 
(20) for the speed u>2- From the next order follows that 5 = 0, that is the perturbations due 
to the coefficients N±, N2, N3 and N± do not modify the speed of the wave while they modify 
the coefficients k„ related to the attenuation in the form 



H _ N 2 (w 2 2 - Vj) + w 2 {pe L N 4 V 2 2 + N lW2 + N 3 v Q - N 4X o) . , 

2wl ■ [ ] 

Anyway, this modification will be small, because w\ — V 2 2 is small and the coefficients Ni are 
also small in the situation considered. 

Now, we consider the case in which n and the initial heat flux qo are orthogonal, and in 
particular we choose n = (0,0, 1) \7\. Through this choice the system (15) becomes 

—pcytjjT — peiuL + Kq 3 = 0, 
(_ w _ iNl )q 3 + K( f + K X0 L = 0, 

(-w - iN 2 )L + Ku q 3 + iN 4 qt = 0, (23) 
l-oj - iN x )qi - iN 3 L = 0, 
(-w - iNi)q 2 = 0. 

If now we make the same assumption as the previous case, i.e. high frequencies uj and high 
wave number k r , the last system becomes 

—pcyioT — pe^L + k r q~3 = 0, 
-ujq 3 + k r ( T + k r xoL = 0, 

-uiL + k r v q 3 = 0, (24) 
-ujq\ = 0, 
-uq2 = 0. 

The dispersion relation of this system is found by setting to zero the discriminant of the matrix 
associated to (24), from which the same speed W2 = oj/k r of the wave of the previous case is 
obtained 

wl = V 2 2 (l-v pe L ) + vl o . (25) 

Now, following the same procedure of the previous case, we consider the quantities Ni, N2, 
N 3 and -/V4 as small perturbations of the physical system in such a way the speed of the wave 
is modified by a quantity 5 and the imaginary part cannot be negligible. Therefore, assuming 
that the new speed of the wave has the form w = co/k r = u>2 + 5, it is found the relation (25) 
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for the speed W2, at the lowest order, and 5 = and 

. _ N 2 {wl-Vi)+N lW l 

at the successive order. Note that also in this case the perturbations do not modify the speed 
of the waves but only the dissipative term k s . From a comparison between the two relations 
of k s , (22) and (26), one may note that the quantities N3 and -/V4 take part only in (22), that 
is, only when the direction of the wave propagation is collinear to the initial heat flux. In 
, ■ , , || , 1 , (p^V 2 2 -Xo)A r 4+A r 3Xo 

particular, one may write k s = kt~ + - 7^2 • 

In summary, we have shown that at high enough frequencies the dynamics of inhomoge- 
neous vortex tangles shows a crossover from a diffusive to propagation behavior. We have 
studied the consequences of these contributions on the speed of high-frequency second sound 
in (20) and (25) and on the attenuation of longitudinal and transverse second sound in (22) 
and (26). The origin of the longitudinal density waves predicted by the equation (2) is to 
be found in the vortex contribution to the total thermodynamics pressure of the system, as 
described in (10). This is different from the transverse elastic waves known as Tkachenko 
waves found in rotating vortex arrays. An evaluation of the speeds shows that the speed of 
the second sound, according to (20) and (25), will not be much influenced by the presence of 
the vortex waves; in contrast, the speed of vortex waves, which would be rather small in the 
absence of second sound, is much increased in the presence of high-frequency second sound, 
because in these circumstances they propagate at a common speed (20) or (25). The results 
for the attenuation coefficient according to (22) and (26) are also interesting, because it turns 
out that the attenuation of second sound at high frequency will be very low. This is in contrast 
with what happens at low frequency, or when the vortex tangle is assumed as perfectly rigid, 
not affected by the second sound, in which case the relative motion of the normal fluid with 
respect to the vortex lines yields an attenuation which allows to determine the vortex line 
density L of the tangle [HE]- However, the wave character of vortex density perturbations at 
high frequency makes that vortex lines and the second sound become two simultaneous waves 
with a low joint dissipation, in the first-order approach. Thus, from the practical point of 
view, it seems that, at high frequency, second sound will not provide much information on 
the vortex tangle because the influence of the average vortex line density Lq is small both in 
the speed as in the attenuation. We have also obtained an explicit expression for the vortex 
diffusion coefficient (13) in terms of k and Buy , which is somewhat analogous to the Einstein 
expression for the diffusion coefficient, but with the quantum n of the turbulence instead of 
the thermal energy and the friction coefficient Buy in the denominator. This may be useful 
for studies on hydrodynamics of vortex tangles [3j IU [7] . 
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